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The recent discovery of the so-called Kuiper belt objects has 
prompted the idea that these objects produce dust grains that 
mas contribute significantly to the interplanetary dust popula- 
tion. In this paper, the orbital evolution of dust grains, of 
diameters 1 to 9 jjl m, that originate in the region of the Kuiper 
belt is studied by means of direct numerical integration. Gravi- 
tational forces of the Sun and planets, solar radiation pressure, 
as well as Poynting-Robertson drag and solar wind drag are 
included. The interactions between charged dust grains and 
solar magnetic field are not considered in the model. Because 
of the effects of drag forces, small dust grains will spiral toward 
the Sun once they are released from their large parent bodies. 
This motion leads dust grains to pass by planets as well as 
encounter numerous mean motion resonances associated with 
planets Our results show that about 80% of the Kuiper belt 
grains are ejected from the Solar System by the giant planets, 
while the remaining 20% of the grains evolve all the way to 
the Sun. Surprisingly, the latter dust grains have small orbital 
eccentricities and inclinations when they cross the orbit of the 
Earth. This makes them behave more like asteroidal than com- 
etary-type dust particles. This also enhances their chances of 
being captured by the Earth and makes them a possible source 
of the collected interplanetary dust particles; in particular, they 
represent a possible source that brings primitive/organic materi- 
als from the outer Solar System to the Earth. 

When collisions with interstellar dust grains are considered, 
however, Kuiper belt dust grains around 9 fim appear likely 
to be collisionally shattered before they can evolve toward the 
inner part of the Solar System. The collision destruction can 
be applied to Kuiper belt grains up to about 50 fim. Therefore, 
Kuiper belt dust grains within this range may not be a signifi- 
cant part of the interplanetary dust complex in the inner Solar 
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1. INTRODUCTION 

Asteroids and comets have long been considered the 
two major sources thaCcontribute dust particles to the 


Solar System interplanetarv space (e.g.. Leinert and Grun 
1990. Levasseur-Regourd c: al. 1991). However, since the 
discovery of the first tram-neptunian object (Jewitt and 
Luu 1992. 1993) and the still ongoing discoveries ot such 
objects (see Jewitt and Luu 1995 for an up-to-date sum- 
mary), it has been proposed that small dust particles pro- 
duced from these so-called "Kuiper belt objects” may con- 
tribute significantly to the interplanetary dust population 
that constitutes the zodiacal cloud (e.g., Flynn 1994). They 
may also be an important source of the interplanetary 
dust particles (IDPs) collected in the stratosphere by high- 
altitude airplanes. The fact that the meteoroid detectors 
on board Pioneer 10 and 1 1 recorded impacts at a nearly 
constant rate up to 18 AU (no data available beyond that) 
from the Sun (Humes 1980) also suggests a possible dust 
grain source in the outer Solar System. 

Three key questions must be answered before we can 
estimate the contribution of Kuiper belt (KB) dust grains 
to the interplanetary meteoritic complex: (1) What is the 
dust production rate in the region of the KB? (2) Can 
those grains evolve all the way to the inner part of the 
Solar System? (3) What is the collisional lifetime of those 
dust grains when interstellar dust particles are considered? 
Question 1 requires the knowledge of, in addition to the 
total population of objects in the KB. for example, the size 
distribution, composition, and strength of the KB objects. 
Based on the study by Jewitt and Luu (1995), there are 
—35,000 objects with diameters larger than 100 km be- 
tween 30 and 50 AU. Recent observations from the Hubble 
Space Telescope (HST) also suggest that there are more 
than 2 X 10 8 Halley-sized objects in the KB (Cochran 
et al. 1995). This population is. in principle, capable of 
producing a large amount of dust (for comparison, there 
are only about 250 asteroids with diameters larger than 
100 km in the main asteroid belt). A recent collision model 
of the KB objects also predicts a total dust production rate 
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from 9 X 10 x to 3 X 10" g sec 1 (Stern 1993. 1996). This 
dust cloud may even he visible in the infrared wavelength 
(Backman et al. 1995: Stern 1996). Compares the dust 
production rate in the KB with the 10' g sec 1 dust needed 
to maintain the zodiacal cloud in a quasi-equilibrium state 
(Whipple 1967). it is obvious that the importance of the 
KB dust grains in the interplanetary dust complex cannot 
be overlooked. Based on that, our objectives in this paper 
are to answer the second and third questions with direct 
numerical integrations of the orbital evolution of dust 
grains and apply a simple collision model based on the 
best available interstellar dust grain spatial density (flux) 
data. In our numerical calculations, we simulate the actual 
orbital evolution of the dust grains by including the gravi- 
tational perturbations due to the Sun. seven planets 
(Mercury and Pluto are not included), radiation pressure 
Poyn ting- Robertson (PR) drag, and solar wind drag' 
We do these calculations for dust grains ranging from 
1 to 9 Mm in diameter. Our collisional model is~ based 
" n lhe recent measurement of the mass and flux of 
mterstellar Just grains by the I 'ly^s spacecraft dust 
detector near the orbit ot Jupiter iGrun et at 1993) 
We assume that this interstellar dust spatial density' is 
the same throughout the Solar System and remains 
constant in time. We neglect the effect of the fluctuating 
solar wind magnetic field force on these interstellar- 
particles even though, for the smaller grains, it may often 
exceed the solar gravitational force. 

Because of PR and solar wind drag, small dust grains 
will spiral toward the Sun once thev are released from 
their large parent bodies (e.g.. Wyatt and Whipple 1950, 
Burns et al . 1979). This orbital evolution leads dust grains 
to pass by planets as well as encounter numerous mean 
motion resonances (MMRs) associated with planets. Based 
on our results, trapping into exterior MMRs with Neptune 
and. to a lesser extent, Uranus dominates the orbital evolu- 
tion of those dust grains. Secular resonances are related 
to the precessions of orbits of the dust and planets involved 
(e.g.. Knezevic et al. 1991). When a small dust grain is 
spiraling toward the Sun, the time it takes to pass by a 
secular resonance location is short compared with the or- 
bital precession time scale of planets. Unless a dust grain 
is trapped in an MMR. the effect of secular resonances is 
not important. Gravitational scattering by Saturn or Jupi- 
ter usually controls their final fate. Approximately 80% of 
the dust grains in our numerical calculations eventually 
were scattered out of the Solar System before reaching the 
inner part of the Solar System. The remaining 20% of the 
dust grains were able to pass bv Saturn and Jupiter and 
completed their journeys all the wav to the Sun 
When collisions with interstellar dust grains are consid- 
ered. however. KB dust grains around 9 /on appear likely 
to be collisionally shattered before they can evolve toward 
the inner part of the Sofar System. Thus, collisions may 


limit the role of KB dust grains in contributing to the 
zodiacal light. 

® ur numerical calculations are described in Section 2. 
The orbital evolution of several typical dust grains is pre- 
sented in Section 3. A statistical description^ of results is 
in Section 4. Our collisional model is in Section 5. The 
significance of KB dust grains in the inner Solar Svstem is 
discussed in Section 6. 


2. NUMERICAL CALCULATIONS 

The orbital evolution of dust grains in the Solar System 
is affected by the radial solar radiation pressure force, bv 
gravitational forces due to the Sun and planets, and by PR 
and solar wind drag forces. The equation of motion of a 
dust grain with geometric cross section A and mass m 
moving under the influence of the gravitational forces of 
the Sun. seven planets, the radiation pressure force, and 
the drag forces can be written as (e.g.. Burns et al. 1979 
Oustafson ! 9M4 ) 


(j\i m (jM,m 
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where v is the velocity vector of the particle. r s is its unit 
position vector with respect to the Sun. and M, is the Sun's 
mass, while M n and r n are the mass of the « th planet and 
the position vector of the particle with respect to that 
planet. 5, c, and Q p! are the solar energy flux density, speed 
of light, and radiation pressure coefficient, respectively. 
The velocity-independent part of the third term at the 
right-hand side of Eq. (1) is the radiation pressure force 
term, whereas the velocity-dependent parts are the drag 
terms. Solar wind drag is caused by the interaction between 
solar wind ions and the dust particle, and the net effect is 
similar to that of PR drag. The ratio of solar wind drag to 
PR drag is represented by sw in Eq. (1). We have assumed 

uin!? tu °' 35 based on ** recent emulation (Gustafson 
1994). The radiation pressure coefficient, Q pr , depends on 

the properties (density, shape, size, etc.) of the particle. 
Traditionally, a dimensionless quantity, fi , is introduced to 
specify the effect of radiation pressure and PR drag. It is 
defined as 


£ = radiation pressure force _ SA Q pr r; 

solar gravitation force ~ cGM.m 

In our numerical models, we study dust grains having B 
equal to 0.4. 0.2. 0.1. and 0.05. The corresponding sizes 
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(diameters) for these four (3 values are 1, 2, 4, and 9 /xm, 
respectively. The conversion from (3 to size is based on the 
assumption that all particles are spherical “astronomical 
silicates'’ (Draine and Lee 1984) with 2.7 g cm' J density 
and using Mie theory to integrate over the solar spectrum 
(Gustafson 1994). Particles smaller than 1 jam usually do 
not remain in bound orbits once they leave their large 
parent bodies, because of radiation pressure. Particles 
larger than 9 /xm usually do remain in bound orbits and 
spiral toward the Sun under drag forces: however, due to 
the long integration times involved (it took approximately 
10 months to complete our numerical simulation on 1- to 
9-/xm particles), we have not modeled the orbital evolution 
of particles larger than 9 gm. Our model serves as a first 
step to study small KB dust grains. The evolution of much 
larger dust grains remains a project for the future. 

We assume dust grains are released (with zero relative 
velocity ) from parent bodies located at 45 AU, with proper 
eccentricities and inclinations equal to 0.1 and 10°. respec- 
tivelv. These are consistent with the observations and sta- 
itsticai studies which show KB objects are in low -eccentric- 
itv oroits with inclinations less than 12 (Jewitt and Luu 
1995. Cochran et ai. 1995). The initial forced elements are 
calculated by the secular perturbation theory (Brouwer 
and Clemence 1961) using the osculating orbital elements 
of seven planets. The initial proper longitudes of ascending 
node, proper longitudes of the pericenter. and mean longi- 
tudes of dust grains (before they left their parent bodies) 
are each randomly chosen between 0° and 360° to simulate 
their collisional origin. The forced and proper orbital ele- 
ments are combined to form the initial osculating orbital 
elements of the parent bodies. When a dust grain leaves 
its parent body, it immediately “sees’' a less massive Sun 
than does the parent body, due to solar radiation pressure 
on the particle. This changes the semimajor axis, eccentric- 
ity, and sometimes the longitude of pericenter of the dust 
grain. All these effects are included in our calculations. 

The interaction between charged dust grains and the 
solar wind-associated magnetic field is not included in our 
model. The uncertainty involved in the structure, as well 
as the time variation, of the field makes it extremely diffi- 
cult to model. At 1 AU, the Lorentz force, F L , on a typical 
l-Lim dust grain, charged to 10 V, is less than 29c of the 
gravitational force due to the Sun, F g , (Leinert and Griln 
1990). For a 10-/um grain, F L is almost four orders of magni- 
tude smaller than F g at 1 AU. The ratio of F L over F g slowly 
increases with increasing heliocentric distance; however, as 
a dust grain moves around the Sun, solar rotation brings 
different sectors of the spiral structure of the magnetic 
field past the dust grain (e.g„ Schwenn 1990, Mariani and 
Neubauer 1990). Together with the fact that the polarity 
of the solar magnetic field changes on a 22-year cycle and 
the charge on a dust grain changes in an unknown fashion 
with time, the net effecUof the solar magnetic field on a 


dust grain’s orbital evolution may not be very significant. 
The actual answer, however, remains an open question. 

The interactions between a charged dust grain and the 
magnetospheres of giant planets are also not included in 
our calculations. This is because ( 1 ) the solar wind pressure 
modifies outer planetary magnetospheres in an unpredict- 
able way over time; (2) planetary dipole fields are not 
aligned with planetary spin axes and the spin rotation 
phases will be totally uncertain in our calculations; (3) like 
the Earth’s magnetic field, planetary fields may change 
with time; and (4) dust grain charge in a planetary magneto- 
sphere is uncertain. Ice grains and their associated evapora- 
tion are not considered in our model; we model only silicate 
dust grains of density 2.7 g cm L 

For a given /3, the equations of motion of seven planets, 
which gravitationally interact with one another, and 20 
dust grains are integrated by the implicit Runge-Kutta 
integrator. RADAU (Everhart 1985), on an IBM ES 9000 
supercomputer and an HP 9000 715/75 workstation. We 
follow the orbital evolution of each dust erain until it either 
spirals into the Sun or is elected Irom the Solar System 
(due to gravitational scattering by planets). 

3. EVOLUTION OF DUST GRAINS 

The evolution in semimajor axis, eccentricity, and incli- 
nation of one of the 9-^rn dust grains is shown in Fig. 1. 
This dust grain started its journey from about 47.5 AU 
(the effect of solar radiation pressure force increases the 
semimajor axis of a dust grain on its release from the large 
parent body). It passed through several exterior MMRs 
with Neptune before it was trapped in the 4 : 3 MMR with 
Neptune (at 35.2 AU) for about 3.6 myr. After it escaped 
this resonance, it was trapped briefly in the l : 1 MMR with 
Neptune (at 29.6 AU). It continued its journey toward the 
Sun after that and eventually was scattered out of the Solar 
System by Jupiter. 

The next example is a 4-/xm dust grain (Fig. 2). This 
dust grain started its journey from about 50 AU. It was 
trapped in the 3:2 MMR with Neptune (at 38.1 AU) for 
about 16 myr. While trapped in this resonance, its eccen- 
tricity increased to about 0.37 before it escaped via a close 
encounter with Neptunfc. "After it escaped the resonance 
via a close encounter with Neptune, it was trapped briefly 
in the 1 : 1 MMR with Uranus (at 18.6 AU). After it escaped 
the last resonance, this particle was able to continue its 
journey toward the Sun until it was removed from the 
Solar System by a close encounter with Jupiter. 

The next example is a 2-/xm dust grain. It started its 
evolution from about 61 AU. due to radiation pressure. 
After passing several MMRs with Neptune, it was trapped 
in the 7:5 MMR with Neptune (at 35 AU) for about 1.8 
myr (Fig. 3). After it escaped that resonance, it was able 
to pass by the four giant planets and made its journey all 
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MMR with Neptune (at 34.6 AU) and then in a 1 : 1 MMR 
with Neptune (at 29.1 AU). Its final resonance trap was 
in a 1:1 MMR with Uranus for about 8.2 myr. In both 
long-duration resonances, the inclination of the dust grain 
tended to decrease with time. Based on previous work 
(Liou et al . 1995b), a secularly decreasing inclination might 
be expected. After it left the last resonance, it continued 
its journey toward the Sun until it was gravitationally scat- 
tered out of the Solar System by Saturn. 

The effect of solar radiation pressure is to change the 
effective mass of the Sun as “seen " by dust grains. This 
changes the semimajor axis, eccentricity, and sometimes 
the longitude of perihelion of a dust grain on its release 
from the large parent body. It also modifies the location 
of a p : q MMR with a planet. The new MMR location is 
given by 


«pi(l - fiy ' 


\P ' 

■Jl. 


= a P ,( 1 - fiV 'K ? '\ 


(3) 


-. -'0 

A ; 

5 10 15 

Time (Myears) 

FIG. 1. Orbital evolution of a 9-/xm KB dust grain. It started its 
journey from 47.5 AU due to the effect of radiation pressure. It was 
trapped in the 4 : 3 MMR with Neptune for about 3.6 myr and was removed 
from the resonance before its eccentricity reached the possible maximum 
value of 0.31. It continued its journey afterward, and was trapped in the 
1 : 1 MMR with Neptune for about 650,000 years. Eventually this grain 
passed by Neptune. Uranus, and Saturn before it was ejected from the 
Solar System by Jupiter. To avoid possible confusion with dust grains 
that spiral into the Sun, we plot the semimajor axis of a grain as positive 
infinity when it becomes hyperbolic. 
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the way to the Sun. When this dust grain crossed the orbit 
of the Earth, its eccentricity was surprisingly small, about 
0.08. Its inclination was about 15°. 

The next example is a dust grain (Fig. 4). The 

radiation pressure force pushes its semimajor axis to about 
130 AU alter its release from the parent body. It was 
trapped in a 2 : 1 MMR with Neptune (at 40.3 AU) for 
about half a million years. After it escaped the resonance, 
it evolved all the way to the Sun. When this dust grain 
passed by the orbit of the Earth, its eccentricity and inclina- 
tion were 0.132 and 16.6°, respectively. 

Another example of evolution is shown in Fig. 5. It is a 
4-/xm dust grain. This is one of the few grains in our samples 
that was trapped in two relatively long resonances. It 
started from about 51 AU and was trapped in a 7 : 4 MMR 
with Neptune (at 42.2 AU) for about 6.1 myr. After it 
escaped the resonance, ft was trapped briefly in a 13: 10 



Time (Myears) 


FIG. 2. Orbital evolution of 4-/u.m KB dust grain. It began its evolu- 
tion from about 50 AU. It was trapped in the 3:2 MMR with Neptune 
for about 16 myr. While trapped, its eccentricity increased to about 
0.37 before being removed from the resonance by close encounters w-ith 
Neptune. This 0.37 is consistent w-ith the maximum eccentricity predicted 
by Eq. (4). It was later briefly trapped in the 1 : 1 MMR with Uranus and 
eventually ejected from the Solar System bv Jupiter. 
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by the planet. A dust grain’s spiral motion toward the Sun 
can be interrupted by getting trapped into one or more 
MMRs with that planet. In a planar circular Sun-planet- 
dust three-bodv system, if the drag rate is small and the 
planet is massive enough such that the evolution o t e 
grain is in the adiabatic regime, the capture probability 
into a given exterior MMR can be estimated (e.g., Henrard 
1982. Beauge and Ferraz-Mello 1994). Otherwise, numeri- 
cal integration of the equations of motion of the system 
has to be performed. As /3 increases, dust grains tend to 
become trapped in exterior resonances closer to the planet 
(Sidlichovsky and Nesvorny 1994). Once the grain is 
trapped in an exterior resonance, its eccentricity increases 
and may eventually reach a maximum value given by 


2 + 3e 

2(1 - e 2 ) 3 


;= K. 


(4) 


This expression gives the value of the eccentricity required 
bv the condition of a quasi-stable trap ( Beauge and Ferra/- 
Mello 1994. Sidlichovsky and Nesvorny l l >94. Liou er ai 
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FIG. 3. Orbital evolution of a 2-jim KB dust grain, it was trapped 
in the 7 5 MMR with Neptune for about 1.8 myr. While trappe . i s 
eccentricity increased to about 0.34 before being removed from the reso- 
nance by close encounters with Neptune. The 0.34 is consistent with 
maximum eccentr.c.ty predicted by Eq. (4). After the grain escaped the 
resonance, it was able to pass by all the planets and completed its journey 
to the Sun. Even with the large variations in eccentricity and inclinati 
over its evolution, this dust grain finally had a small eccentrici y 
and inclination (15°) when it crossed the orbit of the Earth. 


where is the semimajor axis of the planet, and p and q 
are two integers that specify a particular. p:q, resonance. 
The parameter K . defined as the ratio p/q. is farger than 
1 for an exterior MMR. less than 1 for an interior MMR. 

and equal to 1 for a 1 : 1 MMR. 

In a Sun-dust two-body system, the effect of PR drag 
(and similarly solar wind drag) is to take away orbital 
energy and angular momentum of the dust grain. While 
the orbit of the grain circularizes and spirals toward the 
Sun. its inclination (with respect to an arbitrary reference 
plane) remains unchanged (e.g.. Wyatt and Whipple 950. 
Burns et al 1979). When a third body (planet) is involved, 
however, the evolution of the grain becomes more compli- 
cated In addition to the short-period oscillatory as well 
as the secular perturbations from the planet (e.g., Brouwer 
and Clemence 1961), two important dynamical effects are 
brought into the system. They are the MMR perturbations 
by the planet as well as sihgle-event gravitational scattering 
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FIG. 4. Orbital evolution of a l-pm KB dust grain. It started from 
130 AU and was .rapped in a 2:1 MMR with Neptune for half a million 
years. After it escaped the resonance, this gram evo ve . 

the Sun. Its eccentricity and inclination were 0.132 and 16.6 when 
passed by the Earth. 
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FIG. 5. Evolution of a 4- M m KB dust gram that was trapped in two 

!rapped S m a b nC , e h : * ^ W '‘ h NeptUne 3nd a 1:1 with Uranus. While 
pped irt both resonances, the inclination of the grain went throueh 

secu ar. though oscillatory, decreases due to PR and solar wind drag. 

MVfpl : howev f r ’ due to the unstable nature of the exterior 
MMR trap with PR drag (Gomes 1995), dust grains eventu- 
allv eseape the resonance via an ever-increasing amplitude 
of the resonance angle oscillation and. finally, via close 
encounters with the planet. Afterward, they cither con- 
mue their journey toward the Sun or leave the Solar Svs- 

e T MMR Perb b 01 ^ 0rbitS ' If 3 dUSt grain is tra PP ed in a 

1. 1 MMR with the planet, its Jacobi "constant” (when 
ere is no drag this is a constant of the motion in the 

n , t oIrr b ° dy Pr0blem ' e 'S” see Brouwer and 

Clemence 1961) increases with time, which forces the orbit 

to expand from a tadpole to a horseshoe (in the rotating 
reference frame), and. finally, close encounters with thf 
planet eventually remove the dust grain from the resonance 
(L.OU a a. 1 995b). For any dust grain that becomes trapped 
."to an in, error MMR with a planet, trapping is a £ P av , 
unstable and temporary (Liou et at. 1995b) 

In a three-dimensional Sun-planet-dust system, where 
the orbi of the dust grain is inclined with respect to the 
orbital plane of the planet, the dynamics become somewhat 
compheated. The inclination of a dust grain trapped in an 
MMR tends to be damped due to PR and solar wind drag 
'Liou „ at. 1995b); however, secular resonances (when 


in° ‘J h h , P 3net 1S lnvolved ) and gravitational scatter- 
mgs with the planet can also change the orbital inclination 
of the dust gram. They may well dominate drag effects 
3nd inCrease inc 'ination. High inclinations will increase the 
encoomcr velocities of KB dost grains with the Earth 

fob Sy' t^Th' ' huir,0urnev a " lh<! >" 'be inttei 
Solar System. This is an important factor when isidering 

the possible contribution of the KB dust gi ms to the 

collected IDPs. We discuss this in more detail in Section 6 

cated^han " ^ ^ " ine P ' a " etS ’ is more 
cated than a simple three-body system: nevertheless, when 

• 3 d grain 1S lra PP ed in an MMR with a planet the 

dominant perturbation, besides the Sun’s influence, is due 
to that planet A three-body approach to understand the 
namics involved is a good approximation. The evolution 
of dust grains in our numerical simulations are in qualita- 
^ 3 ~ ( WIth . ,heorv ' While trapped in exterior 
‘ , R J’ , h dust grain s eccentricity increases to about the 
predicted maximum value given hv Eq. (4), unless close 
encounters with the trapping planet interrupt the trappmu 
I rocess. or another plane, ,v a .so ,n some resonance w„h 
Jhe grain. I he orbital inclination ot a dust grain in a long 
v "ance can be seen to go through a large short-term 
variation added to a slow secular decrease when a secular 
resonance is not involved. r 

Fip A b U -?h SUai eV ” ,Ution of the KB d “st grain is shown in 
Fig. 6. This is a 9-fim dust grain. Shortly after it started 
its journey from 48 AU, it was trapped in a 7:4 MMR 
with Neptune (at 42.97 AU) for nearly 150 myr. While it 
was trapped, its inclination increased from about 10° to 
more than 40°. This is an example of how a secular reso 
nance can change the inclination of an orbit. At 42 97 A U 
the precession rate of the longitude of the ascending node’ 
ot the dust grain is half the e.genfrequenev of the system 
hat controls the precession rate of the proper longitude of 

sLlv C hm h ng H n ° de ° f NeptUne - Its ecc cntricity increased 
, - 1 had not rea ched the maximum value 0.43 given 

by Eq. (4). The high inclination mav be responsibl! for 
such an unusually long trap. 

4. STATISTICAL DESCRIPTION OF DUST GRAINS 
COMING FROM THE KUIPER BELT 

In our numerical calculations, we calculated the orbital 
evolution of 80 different KB dust grains; there were 20 
grains at each grain diameter of 1. 2. 4, and 9 urn This 
covers approximately three orders of magnitude in terms 
o the particle masses. Surprisingly, their overall evolution 
patterns are very similar. Most of the grains are trapped 

Nenr. 10ns M years ,n ° ne ° f the exterior MMRs with 
Neptune or Uranus. When they escape the resonances 

they are nearly always able to pass by Neptune and Uranus 

and continue their journey toward the Sun. Only three 

dust grains out of the total of 80 samples were ejected 
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(e.g., Opik 1963. Dohnanyi 1978). About three-quarters 
of all sample dust grains were ejected from the Solar Sys- 
tem by either Saturn or Jupiter. Thus, these two major 
planets are primarily responsible for preventing most 
(75%) of the KB dust grains from evolving in toward the 
inner part of the Solar System. 

About 20% of the sample dust grains in each size range 
was able to pass by Saturn and Jupiter. Once dust grains 
passed by Jupiter, they evolved all the way to the Sun. 
None of the dust grains were ejected from the Solar System 
by any terrestrial planets. 

The KB dust grains that do enter the inner Solar System 
(inside Jupiter) have small eccentricities and inclinations 
when they cross the orbit of the Earth. Figure 8 shows 
their eccentricities and inclinations when their semimajor 
axes reach 1 AU. Since PR and solar wind drag reduce 
the eccentricities of dust grains as they spiral from the 
outer Solar Svstem toward the Sun we would expect them 
to have small eccentricities, when no planetarv perturba- 
tions are considered, as they approach the Earth When the 
resonances and eravitationai scatterings that jramattcaiK 
aifect the dust grains' eccentricities and inclinations 
throughout their orbital evolution (Figs. 1-6) are consid- 
ered, however, it is quite unexpected that they approach 
the Earth with asteroidal-type orbits. The low eccentricities 


FIG. 6. Unusual evolution of a 9-p.m KB dust grain. It was trapped 
in the 7:4 MMR with Neptune for nearly 150 myr. At the resonance 
location, the precession rate of the longitude of ascending node of the 
dust grain is two times slower than the proper precession rate of Neptune's 
orbital plane. This secular resonance pumped up the inclination of the 
grain to more than 40°. and overwhelmed drag forces that tend to de- 
crease inclination. 


from the Solar System by Neptune. Trapping intervals due 
to Jupiter and Saturn are short, by comparison. The statis- 
tics are shown in Fig. 7. 

The consequence of dust grains being trapped in exterior 
MMRs with Neptune (or with Uranus) is the formation of 
dust rings like the ring of dust of asteroidal origin near 
the Earth’s orbit (Jackson and Zook 1989, Dermott et at. 
1994b, Reach et al. 1995). Because of their large eccentrici- 
ties, these rings spread out over quite large ranges, approxi- 
mately from 20 to 40 AU. With their large distances from 
the Sun and Earth, it is highly unlikely that these rings are 
telescopically detectable. 

After dust grains pass by Uranus, gravitational scatter- 
ings due to Saturn and Jupiter dominate their final destiny. 
Jupiter and Saturn do not appear to be able to trap dust 
grains in their MMRs in quasi-stable orbits inside 20 AU. 
None of the sample grains were trapped longer than 
500,000 years with either planet. The large masses of Jupi- 
ter and Saturn, however, make their gravitational scatter- 
ings more effective than those of Neptune and Uranus 



FIG. 7. Summary of the evolution of KB dust grains. Even though 
their sizes are different (1 to 9 /urn in diameter. labeled above the right 
histogram), the overall patterns are very similar. Most of the dust grains 
were ejected from the Solar System by Jupiter or Saturn About 20* ot 
the grains were able to evolve all the way to the Sun. 
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Eccentricity 

their eccentricities and incIinattonL n »e y of 'tht asteroidill dust’ "grains" 5 ^ CmCr ' hC S °' ar Sy *‘ em - For sma " »“» <*"»» grains. 


and inclinations of the Earth approaching KB dust grains 
mean they are more likely to be captured by the Earth 
and to survive atmospheric entry than most cometary dust 
grams. After these dust grains passed bv the Earth, they 
continued their journey toward the Sun. Most of them had 
eccentricities smaller than 0.05 when they were within 0.2 
AL lrom the Sun. Their inclinations remained approxi- 
mately the same as those shown in Fig. 8. 

5. COLLISIONS BETWEEN INTERPLANETARY DUST 
GRAINS AND WITH INTERSTELLAR DUST 

In our numerical calculations, mutual collisions between 
dust grains and collisions with interstellar dust particles 
were not included: however, because of the long evolution- 
ary lifetime of the KB dust grains, collisions may destroy 
some dust grains and interrupt their orbital evolution. We 
shall see that collisions may prevent particles larger than 
a certain size from evolving toward the inner Solar System 
even though, from the dynamical point of view, such parti- 
cles are capable of evolving all the wav to the Sun. 

It is difficult to calculate the real collisional destruction 
ol KB grains due to interstellar dust grains. This is due 


primarily to our limited knowledge of interstellar dust (e.g.. 
its spatial distribution and how it varies with time). In the 
following we estimate the role of collisions based on the 

recent Ulvsses measurement and on a few simple assumn- 
lions. * 

In 199_. the spacecraft Ulysses discovered interstellar 
dust grains while passing the vicinity of Jupiter (Griin et 
at. 1993). At 5 AU. the mean flux is 8 X 10' m - sec’ 1 . 
The masses of those particles range from 10' 15 to 5 x l() 
g with a mean mass of 8 x 10 11 g and thev travel in excess 
ol 26 km sec 1 (Griin etjal. 1993, see also Baguhl et al 
1995). According to Eq. (14) in Grim et al. (1985). all but 
the smallest are easily capable of collisionallv shattering 
9- M m KB dust grains. For the purpose of our estimate, we 
assume the mean flux at 5 AU is the same throughout the 
Solar System and that it does not vary with time. From 
the mean tlux. the average time for one collision. t L . for a 
spherical KB grain with diameter d gun to encounter an 
interstellar dust grain, with diameter d, gm. is 


”, , . HEvears. 

d + </,)- 
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If we lake the mean interstellar dust particle mass to be 
8 X 10 13 a, the mean diameter for such spherical particle 
with density 1 g cm 3 is 1.2 /u m. For KB grains with diame- 
ters 1, 2. 4. and 9 M m . the r c are 104. 49. 19. and 4.8 myr. re- 
spectively. 

Mutual collisions between Kuiper belt dust grains can 
also affect their evolution. But the lack of knowledge of 
their size and spatial density distributions in the region of 
KB and the changes of these distributions as different- 
sized grains evolve differently toward the Sun (different 
starting locations, different drag rates, different resonance 
trapping, etc.) all make collision rate calculations very un- 
certain. Here we simply estimate the mutual collision life- 
time of a dust grain in interplanetary space with grains of 
the same size. The average penetration flux for meteoroids 
larger than 8.3 X 10“ 1,1 g measured by Pioneer 10 between 
5 and 18 AU is about 4 X lO - *’ m : sec' ; (Humes 1980). 
The mass threshold corresponds to an 8.4-^.m-diameter 
spherical grain with the same physical conditions as those 
of the dust grains in our calculation. It we assume that this 
constant flux is the same e\en beyond Al' and that the 
meteoroids are ot the same size. 9 /am. then the average 
time for one dust grain to have a collision with another 
one is 


2523 

(9) 2 


10 6 years = 
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The time scale of Eq. (6) is longer than the one given 
by Eq. (5) for a 9-/xm dust grain colliding with an interstel- 
lar grain. When interplanetary dust grains of smaller size 
(with the lower boundary at 1 jam) are included, r c is 
expected to be shortened: however, all KB dust grains are 
in prograde orbits; the impact velocity between them mav 
not be high enough to always result in a catastrophic colli- 
sion. Therefore, mutual collisions between KB dust grains 
appear to be not as important as collisions with interstellar 
dust grains for the 9-/xm dust grains. It is likely true for 
other sizes also. Based on a previous study of the collisional 
balance of the meteoritic complex, which includes the ef- 
fects of PR drag and the observed size distribution of 
micrometeoroids down to 10' 15 g near 1 AU (Grun et al. 
1985), dust grains of 1 to 9-jum have similar collisional 
lifetimes in such a complex. Whether or not these relative 
collision lifetimes remain valid in the outer Solar System 
remains an open question. 

About 90% of asteroidal dust grains smaller than 10 ftm 
in diameter evolve from 3 AU to the Earth without being 
affected by collisions due to the interplanetary dust com- 
plex (Gustafson 1994). Therefore, collisional destruction 
due to the addition of asteroidal dust grains does not seem 
to affect KB dust grains smaller than 10 |xm in the inner 
Solar System. 

The actual orbital evolution lifetimes ol particles ot dit- 


FIG. 9. Dynamical lifetime (from when they were released from their 
parent bodies to the time they were ejected from the Solar System or 
spiraled into the Sun) of the KB dust grains. Several factors determine 
the dynamical lifetime of a KB gram. They are. how large its semimajor 
axis became after it was released from its parent body, its drag rate, 
resonance trappings, and eravitationai scattering by planets. The 9-fxm 
vtrain from Fig. 6, which has a lifetime of 150 myr. is not included here. 


ferent sizes from our numerical integrations are shown in 
Fig. 9. Several factors determine the evolution lifetime of 
a grain. They include the starting semimajor axis of the 
grain (a function of /8 and where the dust gram was released 
from its parent body s orbit), the drag rate, resonance 
trapping, and gravitational scattering by planets. Based on 
our collision estimates' it appears that collisions due to 
interstellar dust grains and mutual collisions have little 
effect on the orbital evolution of 1-. 2-, and 4-^.m KB 
dust grains; however, collisions seem likely to destroy most 
grains around 9 jam before they complete their orbital 
evolution under the PR and solar wind drag. The interstel- 
lar dust grains detected by Ulysses are capable of causing 
catastrophic destruction to dust grains up to about 50 /xm 
in diameter (Grun et al. 1985). For even larger KB dust 
grains the measured interstellar dust grains appear too 
small to cause catastrophic destruction in any one single 
impact, even though “many small projectile impacts ma\ 
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eventually disrupt the target particle although an individual 
impact would not result in complete rupture'' (Grtin et al. 
1985). Also, as the size distribution of interstellar grains 
is not yet well determined, the (lux of large interstellar 
grams is rather uncertain. Whether or not these larger KB 
dust grains can survive collisional destruction and pass by 
giant planets remains an open question. 


6. KUIPER BELT GRAINS AS A SOURCE OF 
ZODIACAL LIGHT AND 1DPS 

Comets and asteroids have long been considered the 
two major sources that deposit dust particles into interplan- 
etary space which, in turn, give rise to the zodiacal light. 
From the in situ spacecraft experiments, it has been shown 
that the main contribution to the zodiacal light is from 
particles that range from 2 ( ! to 200 gin in diameter (Grim 
ft til. 1985). On the other nand. from the Infrared Astro- 
nomical Satellite (IRASi observations at 25 - M m wave- 
length and from dvnamica. modeling, it was shown that l) - 
.um-diameter uum grains - > m comets .mu .iMomids Ivm 
tit the observational shape 1 : the zodiacal light as a function 
ot ecliptic longitudes (Dermott ft til. 1994a. Liou ft id. 
1995a). In regard to the first model, it is not clear whether 
KB dust grains larger than 50 /im can evolve to the inner 
Solar System. In regard to the second model, it appears 
that KB dust grains probably do not contribute significantly 
to the zodiacal light: however, this is based on our collision 
model which assumes a uniform (does not vary with dis- 
tance from the Sun) and stable (does not vary with time) 
interstellar dust grain flux. In reality the spatial density of 
the interstellar dust grains is certain to varv as the Solar 
System moves around our galaxy. It is possible that at a 
time when the Solar System passes through a lower-density 
(lower than what Ulysses detected) interstellar dust cloud. 
KB dust grains larger than 4 gm may be able to evolve 
all the way toward the Sun and be a significant part of the 
zodiacal light. 

Interplanetary dust particles collected by high-altitude 
aircraft from the stratosphere contain several distinctive 
types of primitive materials I e.g.. Bradley ft al. 1988, Clem- 
ett ft al. 1993). These IDPs range from 1 to 50 in 
diameters. They could be responsible for bringing the first 
organic materials to the Earth and contribute to the origin 
of life. Again, asteroids and comets are thought to be the 
two major sources lor these lDPs(e.g.. Bradley etal. 1988); 
however, we argue here that KB dust grains w ith diameters 
around 1 to4/im could be a significant portion of the 
small (between 1 and 4 /ami collected IDPs. 

Short-period comets are thought to come from the re- 
gion of the KB (e.g.. Duncan et al. 1988. Quinn et al. 1990 
Levison and Duncan 1993). They are KB objects that were 
perturbed bv Neptune and thrown into the inner part of 
the Solar System. They fan bring in primitive materials 


from the outer Solar System: however, their mean orbital 
eccentricity and inclination are about 0.5 and 10°. respec- 
tively. Small dust grains released from these short-period 
comets are likely to have eccentricities larger than 0.5, 
if they are still in bound orbits. Because of their high 
eccentricities, their inclinations will be pumped up while 
spiraling toward the Sun (Liou 1993). In general, this leads 
to a higher Earth-approaching encounter velocity (larger 
than 10 km sec ’. Jackson and Zook 1992), although there 
are a few exceptions (Jackson and Zook 1992). There is 
little gravitational focus enhancement for those with high 
Earth-approaching velocities. It is also difficult for them 
to survive the atmospheric entry heating therebv to fall 
down to the surface of the Earth or to be available for 
capture by high-altitude airplanes (Flynn 1989). 

Small dust grains coming from long-period comets may 
not constitute a significant part of the collected IDPs. This 
is because long-period comets have highly eccentric orbits 
with random inclinations. The gain in total enemy bv ..mall 
dust '.trains when thev are released from ionu-perif'u com- 
. is. due to the etfect ot Miiar radiation pressure, can easiiv 
increase their orbital energy to such an extent that tney 
leave the Solar System on hyperbolic orbits shortly after 
leaving their parent bodies. For those that remain in bound 
orbits, the highly inclined orbits (including retrograde or- 
bits) usually result in a very high Earth-approaching ve- 
locitv. 

Large dust grains coming from the main asteroidal belt 
approach the Earth with typical eccentricities less than 0.2 
and inclinations similar to those of their parent bodies 
(with a mean around 10°). This typically results in a low 
(less than 10 km sec' 1 ) Earth-approaching velocity i Jack- 
son and Zook 1992). Small asteroidal dust grains have 
similar eccentricities and inclinations near the Earth (Liou 
and Zook 1995). The gravitational focusing effect enhances 
their flux near the Earth (e.g., Dohnanyi 1978). A low 
Earth-approaching velocity also makes it easier for them 
than for highly eccentric orbit grains to survive the atmo- 
spheric entry heating. They are certainly a source of the 
collected IDPs: however, they only bring in materials from 
as far as about 4 AU from the Sun. Whether primitive/ 
organic materials at that distance could survive the early 
Solar System T-Tauri w T ind and heating from the Sun is 
not clear. 

When one compares the KB dust grains with either com- 
etary or asteroidal dust grains, it becomes obvious that KB 
grains have certain advantages to make up the carbona- 
ceous unaltered small IDPs on Earth. They approach the 
Earth with eccentricities and inclinations similar to those 
ol the asteroidal dust grains. The consequent low Earth- 
approaching velocities not only enhance their flux but tavor 
their chances ol surviving the atmospheric entry heating. 

I hey are likely to bring in primitive/organic materials from 
the outer Solar System, where primitive/organic materials 
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should be better preserved simply because of their greater 
distance from the Sun. They could be a significant fraction 
of the small (less than 5-^.m-diameter) collected IDPs. 
They could even be a source of the earliest organic material 
on the Earth that gave rise to life. 

Is it possible to distinguish KB dust grains from other 
dust grains in the collected IDPs? If so, we can analyze 
them and learn more about material composition in the 
KB and increase our knowledge of the early Solar System. 
From a dynamical point of view. KB dust grains that cross 
the Earth’s orbit appear just like asteroidal dust grains, 
with similar eccentricities and inclinations when they cross 
the orbit of the Earth. This implies that the atmospheric 
entry heating is similar for both types of grains. The major 
difference between KB and asteroid dust grains is their 
space exposure time. The solar flare track density measure- 
ment (e.g.. Bradley et al. 1984. Sandford 1986) may show 
a difference between them. KB dust grains should have a 
much higher track density than asteroidal dust grains. This 
suggests that among small collected IDPs. those with the 
hianest soiar tlare track densitv mav come trom the KB. 

7. CONCLUSION 

We have completed a study of the orbital evolution of 
micron-sized dust grains (1 to 9 /am in diameter) coming 
from the region of the KB. Our model includes the gravita- 
tional forces of the Sun and planets, as well as solar radia- 
tion pressure and PR and solar wind drag on a dust grain. 
Interactions with solar magnetic fields and the magneto- 
spheres of giant planets are not included in the calculation. 
The results show that a significant portion (about 20%) of 
these dust grains are capable of evolving toward the inner 
part of the Solar System. With the continuous discoveries 
of large objects in the KB, the outer Solar System (from 
30 to 50 AU and beyond) has become a place that could 
potentially supply a large amount of dust grains to inter- 
planetary space. Collisional destruction due to interstellar 
dust grains, however, may limit the significance of KB dust 
grains in contributing to the zodiacal light. Small KB dust 
grains (less than 9 gm) can evolve toward the inner Solar 
System without being seriously affected by collisions. For 
much larger KB dust grains (larger than about 50 /Am), it 
remains a possibility that they may also enter the inner 
Solar System. To understand their exact contribution to 
the zodiacal light and in the collected IDPs. and to possibly 
distinguish them from other dust grains in the collected 
IDPs. more w-ork is needed on the dust production rate in 
the KB. the chemistry in the outer Solar System, as w'ell as 
numerical modeling of evolution of larger KB dust grains. 
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